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ABSTRACT 

There is no accepted mechanism that explains the equilibrium structures that 
form in collisionless cosmological N-body simulations. Recent work has identified 
nonextensive thermodynamics as an innovative approach to the problem. The 
distribution function that results from adopting this framework has the same 
form as for polytropes, but the polytropic index is now related to the degree 
of nonextensiveness. In particular, the nonextensive approach can mimic the 
equilibrium structure of dark matter density profiles found in simulations. We 
extend the investigation of this approach to the velocity structures expected from 
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nonextensive thermodynamics. We find that the nonextensive and simulated N- 
body t' rms distributions do not match one another. The nonextensive v rms profile 
is either monotonically decreasing or displays little radial variation, each of which 
disagrees with the t> rms distributions seen in simulations. We conclude that the 
currently discussed nonextensive models require further modifications in order to 
corroborate dark matter halo simulations. 

Subject headings: dark matter — galaxies:structure — galaxies:kinematics and 



Simulations of structure formation have become well-refined and sophisticated over the 
past few decades. Typically, N-body simulations follow collisionless particles that repre- 
sent dark matter and hierarchically form gravitationally bound structures. These structures 
have several "universal" characteristics, among them their self-similar density profiles (e.g., 
Navarro et al. 1996; Moore et al. 1998; Power et al. 2003). However, as of yet there is no 
accepted physical mechanism for explaining these shared properties, although there is no 
lack of hypotheses (e.g., Lokas 2000; Nusser 2001; Barnes et al. 2005; Lu et al. 2006). This 
work further examines the recently advanced suggestion that nonextensive thermodynamics 
can be used to determine the equilibria (Hansen et al. 2005; Leubner 2005; Hansen & Moore 
2006; Kronberger, Leubner, van Kampen 2006). For simplicity, we deal only with spherically 
symmetric systems in this work. 



The lack of short-range forces in self-gravitating, collisionless systems throws the adop- 
tion of standard thermodynamics into question. Instead, Tsallis (1988) has proposed a 
generalized themodynamical approach for such systems based on the nonextensiveness of 
entropy; i.e., the entropy of a composite system is not simply the sum of the subsystems' 
entropies. In particular, the entropy of a composite system is given by, 



dynamics 



1. 



Introduction 



1.1. Nonextensive Thermodynamics 



Sa+b — Sa + Sb + (1 — q)SASB, 




where S is the entropy, A and B refer to subsystems, and q is a parameter that describes 
the degree of nonextensiveness. When q = 1 the situation is extensive and corresponds 
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to standard thermodynamics. Nonextensive thermodynamics (sometimes referred to as q- 
statistics) has been designed to describe systems with long dynamical memories, as in self- 
gravitating collisionless objects. Although other applications have been discussed by Tsallis 
(1999), specific applications of nonextensive thermodynamics to astrophysical situations have 
been discussed by Plastino & Plastino (1993); Aly (1993); Hansen et al. (2005); Leubner 
(2005), and most recently by Kronberger, Leubner, van Kampen (2006, hereafter KLK). 
One of the more interesting results of the nonextensive approach is that the associated 
distribution function has the same form as for a polytrope (Plastino & Plastino 1993). In 
this case, the analogue to the polytropic index is related to the q value (see § 2.1). Hansen & 
Moore (2006) appeal to nonextensive thermodynamics to interpret the correlation between 
density profile slope and velocity anisotropy observed in N-body simulations. KLK have 
found that the density profiles corresponding to nonextensive distribution functions can be 
found to match those of virialized dark matter halos formed in cosmological simulations. 
These studies have demonstrated that nonextensive thermodynamics is an attractive option 
for theoretically grounding the behavior seen in computer simulations. After a brief synopsis 
of relevant halo properties, we will discuss our findings that nonextensive thermodynamics, 
as it currently stands, does not provide such a basis. 



1.2. Overview of Simulated Dark Matter Halo Properties 

The virialized structures (halos) that are formed in cosmological N-body simulations 
have several properties that appear to be "universal". Density profiles have nonpower law 
shapes where a scalelength divides a steeply declining outer region from a less steeply de- 
clining inner cusp. Most previous works investigating halo formation focus on explaining 
this kind of density profile. However, the velocity structures of the halos are also impor- 
tant. Like the density profiles, N-body halo velocity dispersion distributions v ims (r) have 
nonpower law shapes which peak near the density scalelength. The relevance of be 
illustrated by the quantity p/vf ms which is reminiscent of and shares the dimensionality of 
the phase space density. As first noted in Taylor & Navarro (2001), this quantity follows a 
power law distribution with radius p/v^ ms oc r~ a , where a is a constant that depends on the 
particular density profile. This behavior appears to be linked to the physics of collisionless 
collapse, given that the same power laws result from independent semi- analytical halo for- 
mation models, as shown by Austin et al. (2005). The velocity structures of dark matter 
halos therefore provide additional leverage for understanding the physical mechanisms that 
drive equilibrium halo formation. 
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2. Nonextensive Isotropic Equilibria 

KLK demonstrate that the equilibrium density profiles of sufficiently nonextensive sys- 
tems with velocity isotropy describe the results of N-body simulations well (see Figure 2 of 
KLK). Their work maintains a constant degree of nonextensiveness q. To maintain continu- 
ity with KLK, we will adopt their variable k — 1/(1 — q) as the measure of nonextensiveness. 
We first investigate the predicted velocity dispersion for these models, and then we proceed 
to allow n(r) to vary. 



2.1. Constant k 



Starting with the nonextensive definition of entropy (Tsallis 1988), 



S K 



1-1/ K 1 

Pi - 1 



(2) 



where ks is Boltzmann's constant and the sum is over all accessible states (each with prob- 
ability p^, one can extremize S K under the constraints of constant total mass and energy 
to find the distribution function, as in Plastino & Plastino (1993). Silva, Plastino, & Lima 
(1998) find that the same nonextensive velocity distribution functions can be found using a 
variation on Maxwell's derivation of the standard distribution functions, independent of any 
assumptions about energy. As Silva, Plastino, & Lima (1998) make clear, the nonextensive 
approach demands velocity isotropy. 

Assuming a spherical system and adopting the Poisson equation as the link between 
density and potential, the nonextensive distribution function is obtained (Plastino & Plastino 
1993; Kronberger, Leubner, van Kampen 2006), 

E, 



r(E r 



B- 



1 + 



K<7 



(3) 



where E r = v 2 /2 — ^ is the relative specific energy, v 2 /2 is the specific kinetic energy, 
\1/ = — $ + $ (^ is the specific potential energy and $o is a constant), and a is an en- 
ergy normalization constant. The '±' denotes that k can be positive or negative and the 
normalization constant B reflects that choice. Again, note that this distribution function is 
exactly that for a polytrope (Plastino & Plastino 1993). The nonextensive thermodynamic 
framework provides a physical meaning to the polytropic index; it is a measure of nonex- 
tensiveness. With this distribution function [which is isotropic since / = f(E r )], it can be 
shown that the density is, 
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and the rms velocity distribution is given by, 



,± 



KG 2 



1/2 
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(5) 



where w rms =< v 2 > 1 ^ 2 . We note that Equation 5 is reminiscent of the expression for the 
outer halo velocity dispersion given in Equation 3 of Hoeft, Miicket, & Gottlober (2004). 
It is straightforward to show that these functions satisfy the Jeans equation for mechanical 
equilibrium. Equation 4 can be solved for \1/ and combined with the Poisson equation, 

V 2 * = -AirGp (6) 

to form a second-order differential equation for p that we solve numerically. 

We turn our attention here to the velocity profiles given by Equation 4. We will specif- 
ically consider profiles with k < 0, like those discussed in KLK, for reasons discussed at 
the end of § 2.2. For large values of | system becomes increasingly extensive), 

v Tms approaches a constant as the density profile becomes isothermal, as expected. For more 
modest values of k < 0, v rms is always a continuously decreasing function of radius. However, 
this is not what is seen in the results of cosmological N-body simulations, for which v rms has 
a well defined peak near the scalelength of the halo density profile tq. 

Another view of this discrepancy is given by the phase-space density-like quantity p/v^ ms . 
The solid red lines in Figure 1 represent nonextensive density (panel a), velocity dispersion 
(panel b), and p/v^ ms profiles (panel c) using the KLK values k = —15, cr = 0.12. Note that 
the density and p/v^ ms distributions have been scaled by (r/r )~ 2 to highlight differences 
from this power law profile. Figure 1 also shows p, f rms , and p/v^ ms profiles for isotropic 
Navarro et al. (1996, hereafter NFW) (dashed blue lines) and Navarro et al. (2004, hereafter 
N04) (dash-dotted green lines) models for comparison. In the end, although a nonextensive, 
constant k density profile can be found to mimic the density profiles found in cosmological 
simulations, the corresponding velocity structures do not match. 



2.2. Variable n 

We now relax the constraint of constant k and allow for k(t) distributions. As before, 
Jeans equilibrium is maintained. Varying k introduces derivatives of k to Equation 6, in 
addition to the derivatives of p, allowing us to choose either: A) a p distribution and solve for 
n(r), or B) a k distribution and solve for p(r). We have tested our numerical implementations 
of these approaches for self-consistency. 
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Accounting for radial variations of k is especially pertinent given that it has recently 
been found that there is a direct connection between k (or q) and the slope of a density 
profile (Hansen et al. 2005). The radial changes in density slope present in dark matter halo 
models (e.g., density slopes change from -1 to -3 as radius increases in NFW profiles) then 
demand radial variation of k. 

Following track A dicussed above, both NFW and N04 density profiles lead to derived 
k distributions that are always negative. While we discuss the N04 profile in detail here, 
the NFW results are essentially the same. Figure 2 shows the scaled N04 density profile 
(panel a), three k profiles (panel b), and the resulting v ims distributions (panel c). The solid 
blue, dashed green, and dash-dotted red lines represent solutions when K ini t = —5, —15, and 
—20, respectively. We could not find solutions with positive k values over any section of 
the distribution. The resulting k distributions have roughly sinusoidal shapes superimposed 
on linear trends which vary with the adopted initial k value; for K\ n \t > —17 the linear 
trends have negative slopes, for < — 17 the trends are positive. While these profiles are 
interesting in their own right, the relevance to the current work is that the resulting velocity 
dispersion profiles in Figure 2c remain monotonically decreasing functions, like the constant 
k case (solid line in Figure lb). We conclude that the nonextensive distribution given in 
Equation 3 does not reproduce all aspects of N-body simulation results. 

As track A fails to produce velocity dispersions similar to those in simulations, we now 
discuss a k distribution designed to remedy this problem. We seek a function that is positive 
(negative) for radii smaller (larger) than r . This should give us a dispersion profile that 
peaks near ro- Since k = leads to infinite entropy, we cannot simply transition from positive 
to negative values. These arguments have led us to choose the following form, 

n(r) = —A cothfs log (r/r )], (7) 

where k — ► A as logr/Yo — > — oo, k — > —A as logr/Yo — > oo, and s determines the rate of 
change between A and the infinite value at logr/ro = 0. As many functions can be created 
with similar properties, we discuss this specific form only as an illustrative example. 

The k profile described by Equation 7 is shown in Figure 3a. Where k < 0, we plot 
— log | /c | - A density profile that reasonably approximates those from simulations has the 
following parameter values; A = 25, s = 0.5, and a = 0.12 (Figure 3b). As in Figure 1, 
the solid red, dashed blue, and dash-dotted green lines represent solution, NFW, and N04 
profiles, respectively. The v rras distribution now has the correct qualitative behavior; it is a 
rising function near the center and decreases with increasing r > r (Figure 3c). However, 
it does not quantitatively match simulation results. Specifically, the v rms profile is too flat. 
This difference is also apparent when looking at the p/v^ ms distributions in Figure 3d. The 
P/ V rms profile from the nonextensive approach closely resembles the density, which has a 
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decidedly nonpower law shape. There simply is not enough variation in the w rms distribution 
to compensate the changes in p and produce a power law. Such discrepancies may reflect the 
insufficiency of Equation 7; further work may suggest more appropriate n(r) expressions. 

We note that this use of k > for logr/ro < contradicts the arguments made in KLK 
and Leubner (2005). Leubner (2005) points out that negative k means: 1) that the entropy 
of a composite system is less than in the extensive case, and 2) the heat capacity is negative. 
k > implies a positive heat capacity, indicitive of a system in which self-gravity is not 
important. We do not speculate further on the implications of this interpretation. Positive 
k has been adopted simply to get the v rms profile to have positive slope at small r like the 
simulation results. 



3. Summary & Conclusions 

It has been recently shown that nonextensive thermodynamics can be used to explain the 
density profiles of collisionless dark matter halos formed in cosmological N-body simulations 
(Kronberger, Leubner, van Kampen 2006, KLK). While this is an important step forward 
in understanding the physics governing collisionless equilibria, N-body simulations also find 
universal links between density and velocity distributions, such as the power law behavior of 
P/ V rms i e -9-! Taylor & Navarro 2001). Can nonextensive thermodynamics correctly predict 
the velocity behavior of systems formed in cosmological simulations? 

We find that a constant nonextensiveness parameter k (as in KLK) predicts a veloc- 
ity dispersion profile that is a continuously decreasing function of radius, in conflict with 
simulated profiles. We have extended the KLK approach to include radially varying k dis- 
tributions designed to more closely match the velocity behavior of simulation results. While 
this approach produces qualitative matches between the predictions and simulation results, 
there is no quantitative agreement. We describe four possible solutions to this problem. 

1) Equations 1 and 2 correctly describe long-range collisionless gravitational interactions, 
but cosmological N-body simulations are not purely collisionless close to centers of halos due 
to numerical effects. In this case, a velocity dispersion that increases with radius (for r < ro) 
may indicate that artificial processes that mimic short-range interactions might be at play in 
the simulations. We do not speculate on what these might be (but see El-Zant 2005) or why 
they would lead to velocity dispersions increasing with radius. However, halos calculated 
from semi-analytical collapse models can reproduce both the density and velocity profiles of 
N-body halos, arguing against numerical effects. 

2) The nonextensive approach of Equations 1 and 2 correctly predicts the global (asymp- 
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totic) equilibrium state, but N-body simulations (and possibly, real dark matter halos) set- 
tle to nonlocal or quasi-equilibria through instabilities, such as the radial orbit instability 
(Barnes et al. 2005). 

3) Equations 1 and 2 do not correctly describe collisionless systems with long-range 
interactions, so a different thermodynamic approach needs to be found to describe dark 
matter halo simulations. 

4) The approach taken in KLK and here does not account for velocity anisotropy which 
is present in halos formed in cosmological N-body simulations (Cole & Lacey 1996; Huss et 
al. 1999; Barnes et al. 2005; Hansen & Moore 2006; Hansen & Stadel 2006). Preliminary 
work by the authors that incorporates anisotropy into the nonextensive framework suggests 
that this is not a panacea. This adds to the evidence presented in Hansen et al. (2006) that 
suggests the presence of anisotropy in N-body halos requires that the distribution function 
have a more complex form than that given by Equation 3, pointing strongly towards the idea 
that a different thermodynamic approach must be taken for these systems (point 3 above). 

Regardless of the validity of these points, further investigation of the relationship be- 
tween nonextensive thermodynamics and self-gravitating systems is vital to providing a 
physical understanding, and independent corroboration, of the results of dark matter halo 
formation simulations. 
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Fig. 1. — The constant k nonextensive solutions for density (panel a), velocity dispersion 
(panel b), and p/uf ms (panel c) are given by the solid red lines. For comparison, correspond- 
ing curves for the NFW (dashed blue lines) and N04 (dash-dotted green lines) models are 
also shown. The density and p/t'j? ms profiles have been scaled by (r/ro) 2 to highlight varia- 
tions from this power law distribution. The nonextensive v vms distribution is monotonically 
decreasing, in stark contrast to the behavior of the simulation-based NFW and N04 models. 
This discrepancy is also evident when comparing the decidedly nonpower law shape of the 
solution p/t>r ms profile to the very nearly power law empirical p/ff ms profiles. Note that the 
discrepancy is evident at radii that are well-resolved by simulations. 



-12- 




Fig. 2. — When k is not constant, a k profile can be found by solving Equation 6 with an 
assumed density distribution. The results shown here use the N04 density profile; adopting an 
NFW profile leads to very similar results. Panel a shows the scaled N04 density distribution 
that is shared by all of the solutions. The solution k profiles (with = 0) are shown 
in panel b. Here, the solid blue, dashed green, and dash-dotted red lines indicate solutions 
with K init = —5, —15, and —25, respectively. The velocity dispersion profiles (panel c) are 
montonically decreasing, as in the constant k case. Similarly, these v Tms profiles do not match 
those from simulations. 
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Fig. 3. — For the k distribution described by Equation 7 (panel a), we show the nonextensive 
solutions (solid red lines) for density (panel b), velocity dispersion (panel c), and p/v^ ms 
(panel d). For comparison, corresponding curves for the NFW (dashed blue lines) and N04 
(dash-dotted green lines) models are also shown. As in Figure 1, the densities and p/v^ ms 
profiles have been scaled by (r/Yo) 2 to highlight variations. The solution curve is a decent 
match to the empirical density profiles in panel b. Close inspection of panel c reveals that 
the solution v ims profile does first increase before decreasing, in qualitative agreement with 
the NFW and N04 curves. However, it is quantitatively very different from those curves, a 
fact also reflected in the discrepancies apparent in the p/v^ ms profiles in panel d. 



